The lower bound of barrier-energy in spin glasses: a calculation of the exponent on 

hierarchical lattice 
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We argue that the lower bound to the barrier energy to flip an up/down spin domain embedded 
in a down/up spin environment for Ising spin glass is independent of the size of the system. The 
argument shows the existence of at least one dynamical way through which it is possible to bypass 
local maxima in the phase space. For an arbitrary case where one flips any cluster of spin of size 
I, we have numerically calculated a lower bound to the exponent ip characterizing the barrier one 
has to overcome. In this case ip corresponding to the lower bound calculated on hierarchical lattice 
comes out to be equal to 6 the exponent characterizing the domain wall energy in ground state. 

PACS numbers: 



The slow dynamics of spin glasses [SG] after a quench 
from a temperature greater then T c to one below T c (T c - 
spin glass transition temperature) is to a good extent ex- 
plained under the assumption of droplet theory 0, of 
coherent excitations of spin domains. The aging of a spin 
glass system (below T c ) is understood as a slow increase 
in coherence length Lxit) which occurs by accumulating 
smaller domains into bigger ones and in that way reduc- 
ing the total domain wall length in the system. This, 
very slow non-equilibrium behavior which leads the sys- 
tem to-wards equilibrium is in excellent agreement with 
the scaling law 
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where t is time with tq(T) being the time scale, T is 
the temperature of the system whose energy has been 
calculated in the unit of Boltzmann constant ks and ip 
is the exponent associated with the free energy barrier 
one has to overcome in order to go from one minimum 
to another. It is taken that the barrier energy scales as 
& where I is the system size. The barrier exponent is 
considered to satisfy the relation 9 < ip < d — 1. The 
upper bound comes because of the fact that any com- 
pact cluster of length scale I and energy less than l^ 1 
can be created with a maximum domain wall size 
by simply moving the domain wall through the cluster, 
and the lower bound corresponds to the fact that one 
cannot flip a domain of size I without encountering one 
of size at least 1/2 whose free energy scales as (l/2) e 0. 
The exponent 9 being much smaller than (d — I) char- 
acterizes the free energy of a spin glass domain of linear 
size I. A lot of experimental end a few numerical work 

HSHIlSESllS have shown the value of V> to fall 
within the above mentioned range. The fact that some 
experimentally evaluated value of ip in two and three di- 
mensions is close to its upper limit is recently raising 
some question about the validity of the lower limit of ip 



set by Fisher and Huse Q . In a recent paper it has been 
argued in favor of a lower limit as ip = d — 1 rather 
than ip = 9. An efficient numerical scheme has been pro- 
posed in Ref.-^lJ following some underestimations rele- 
vant to the lower bound to barrier energy. In the present 
work we will mainly follow this numerical procedure (ap- 
plicable for very high dimensions and length scales) but 
to definitely get different results than that in [llj . 

Getting an idea about the lower bound of ip in high 
length scales is a notoriously complicated affair in view of 
the fact that one ideally requires local optimizations of all 
possible spin flipping moves which is an NP-complete [l^ 
process. Let us first present a dynamical way of sequen- 
tial small excitation and subsequent stabilization which is 
able to slowly shrink the domain wall covering an entirely 
up/down spin domain embedded in a down/up spin en- 
vironment. Fro simplicity, consider a rectangle separated 
by a straight line down the middle into two parts. In one 
part there are spins all in up-wards direction and in the 
other part there are spins all in the down wards direction. 
The middle line will be pushed to-wards one edge in such 
a way so that in the end the whole rectangular region is 
occupied by down spins only. In reality the domain wall 
is definitely not a straight line but much thicker than 
but our result can easily be generalized to that. Let us 
think about any sequence of single spin flip sequence that 
pushes the middle st line near its one end to-wards the 
up-spin domain living behind a very small extension of 
the down spin domain. Our changed domain wall i.e. the 
small pushed in part of the straight line, now has many 
spins frustrated on it immediately following the change in 
shape and we stabilize them to attain the ground state in 
the changed configuration of the spin domains. Since, in 
the above mentioned process we have only excited a very 
small portion of the domain wall and moreover, after the 
stabilization in the new state the domain wall length has 
changed only by a very small amount, the barrier energy 
encountered is very small. Now, we repeat the process of 
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pushing the straight line (domain wall) in the same way 
just in the adjacent part of the previously pushed one 
and do the whole process. In the process, still we keep 
having the same stabilized length of the domain wall as 
in the previous state and had only excited a small part 
of it. So in this process we can make the down spin state 
encroach the up spin region totally but never going at a 
very high energy state compared to the initial one. 

Now, to generalize the above mentioned dynamics, 
consider an up/down spin domain completely embedded 
in a down/up spin one having a domain wall in the mid- 
dle which can have a complex shape. Flipping the middle 
up/down spin domain can be achieved by slowly shrink- 
ing the domain wall in the above mentioned way to even- 
tually make it vanish. In such a process since we are 
shrinking the boundary wall of a simply connected do- 
main, we might at worst need to increase the domain 
wall by a very small amount. Moreover, since the ex- 
citations at each step of movement are kept very small 
the barrier state energy should never be very large com- 
pared to the initial state. If the up/down spin domain 
is multiply connected - meaning that there are droplets 
of oppositely oriented spins in the middle - the trick will 
be to flip all the smaller droplets inside following the 
same procedure to make the inner domain walls vanish 
first and then make the bigger domain wall go to zero 
as mentioned above. Thus, we can argue that there is 
at least one dynamical way corresponding to that men- 
tioned above, where the barrier energy should not scale 
as the system size. For such a situation, there is defi- 
nitely a way to bypass the local maximum in the phase 
space when one goes from one region to the other and 
the path corresponds to the lowest barrier path. 

But, definitely the above mentioned way is a too 
stricter one to follow considering thermal motion to pro- 
duce it. In the following we would consider the way that 
involves a distribution of the excitations where larger do- 
mains can be flipped coherently at a time. Here one 
should think about optimizing on the flipping sequence 
of intermediate sized clusters to introduce or abolish a 
domain wall at a length scale /. In the rest of this paper 
we are going to show that the lower bound to ip is equal to 
6 for an Ising spin glass on a Migdal-Kadanof [MK] Hier- 
archical lattice. In what follows we apply renormalization 
group transformations to directly calculate ip at its lower 
limit at various length scales in many dimensions and 
will show the agreement in them at lowest and highest 
length scales (in our numerics) in all dimensions starting 
from d — 2. Keeping in mind that it is an impossible 
task to probe all the relevant single spin flip sequences to 
make sure that one of it passes through the true barrier 
state one has to apply the optimized scheme of sequence 
of cluster flips with a cluster length-scale distribution ap- 
propriate to the given length scale of the system. Let us 
think about optimizing on a sequence of flipping cluster 
of spins before changing the boundary condition which 



will either introduce or remove a domain wall inside the 
system of size I. In such a process we would typically 
encounter the barrier state just before or after we change 
the boundary condition. When we flip a single cluster of 
spins which was previously in ground state in the process 
of flipping a sequence of them before changing the bound- 
ary conditions we suddenly frustrate a lot of spins on the 
previous domain wall of that cluster. In principle, one 
could find a single spin flip sequence for the spins on the 
domain wall before flipping the cluster as a whole which 
would at best lead one to a barrier state corresponding 
to the flip of that cluster with half as high energy. This is 
so because in this case also one would encounter the bar- 
rier state just before or after flipping the cluster followed 
by adjustment of spins on the domain wall. The error 
in this estimate will be of the order of energy associated 
with a single spin flip. Thus, the scope of adjustment 
of spins before flipping an intermediate sized cluster in 
a sequence of cluster flips to get to the barrier state as- 
sociated with change in boundary conditions of a system 
of size I could result in a barrier state half as high and 
will not essentially alter the scaling law. Since the condi- 
tions remain unchanged for spins lying inside the cluster 
(not on or adjacent to domain wall) undergoing flip we 
do not bother about them. Here, we actually consider 
that within cluster adjustments of spins after its flip as a 
whole are much low energy affairs to change our result. 

Our numerical calculations are at T = on Ising spin 
glass. An Ising spin glass has the Hamiltonian 

H = -y]jjjSjSj (2) 

(«> 

and undergoes spin glass transition at nonzero T c in di- 
mensions d > 2. In Eq.2 the sum is over all nearest 
neighbour pairs and the spins S% and Sj can take up val- 
ues ±1. There are good agreements of calculated value of 
9 on M if-hierarchical lattice and square and cubic lat- 
tice in two and three dimensions 0, Q, 0, 0] and a 
hierarchical lattice is widely in use for the numerical in- 
vestigation of spin glasses. The Mi^-hierarchical lattice 
starts by forming in its first level a unit with 2 d ~ 1 parallel 
bridges connecting two end spins which set the bound- 
ary condition. Every bridge is a series connection of two 
bonds with a spin in the middle. In the present case at 
first level we select the bonds randomly from a pool of 
10000 bonds generated from a Gaussian distribution of 
unit width and zero mean. In the next higher level of 
length scales one replaces the previous bonds with ones 
obtained from a renormalized bond distribution which is 
equivalent to replacing each bond in the previous level 
by the whole unit. Thus, in the 1th level of iterations 
the renormalized lattice corresponds to a length scale 2 1 
where each spin in the middle is a domain of intermedi- 
ate size. On such a lattice at T = 0, the renormalization 
group transformation of effective bonds involves getting 
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a contribution from each bridge as 

Ji = sign(JiJ2)min[\Ji\, \ Ji 



(3) 



where J\ and Ji are the bonds connected in series in the 
bridge. The renormalized bond is obtained by summing 
over contributions obtained from all parallel bridges. 

The basic theory we follow in optimizing the middle 
spin flip sequence before we flip the right hand corner 
spin to change the boundary condition is the same as 
in Ref.|ll| apart from the fact that we are using renor- 
malized bonds and are not concentrating to L d_1 nearest 
neighbour spin to the right hand corner spin. Let us ex- 
plain the steps of spin flip sequence following Ref.jllj. 
We would always go from the higher energy minimum 
to the lower one corresponding to the states with paral- 
lel and anti parallel end (corner) spins and will find the 
barrier height from the reference energy of the higher ly- 
ing ground state. Since we are after finding the lower 
bound to the energy barrier the above mentioned move 
is justified. Now consider one of the 2 d ~ 1 bridges that 
connect the end spins. Let us take e 1 as the energy of 
the bridge measured from the reference of the previous 
ground state energy of it when the middle spin is flipped 
first before flipping the right hand corner one. Let e 2 be 
the energy of it when the corner spin if flipped first and 
the middle spin remains as it was and e 3 be the energy of 
the bridge when the middle spin is flipped first and then 
the corner spin is flipped. Taking the number of middle 
spins flipped before we flip the right hand corner spin to 
change the boundary conditions as n, the energy of the 
system just before and after flipping the corner spin be 
respectively 
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and 



E <* 2 + £ 



ei 



(4) 
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k=2 d 



1=1 



Since the energy barrier will be encountered just before 
or after flipping the end spin the barrier energy is given 
by 



E = Max[E\E 2 



(6) 



As we are searching for the lowest barrier energy, we can 
make some underestimations. For a parameter a satisfy- 
ing < a < 1 we can write 

E = Min n [Max(E l ,E 2 )] 
= Min n \aE x + (1 - a)E 2 ] 



= Min n 



5^ (aej 1 + [1 - a] 



ei 



2 d-i 

E <* 2 

k=2 d - 1 -n 



2 a-i 

E 

i=l 



Minliae 1 + [1 - a]e 3 ), (1 - a)e 2 



(7) 



In the above expression setting a — 0.5 means giving 
equal weight to E 1 and E 2 which will actually make the 
barrier energy negative in 2-dimensions This is be- 
cause of the fact that we are always going from higher en- 
ergy minimum to a lower energy minimum and it makes 
e 2 negative more often than not in 2-dimensions. In fig. I 
we have plotted the log of barrier energy E against log 
of system size (I = 2 1 ) for a = 0.9 . The plot shows 
graphs corresponding to dimensions from d = 2 to d = 9. 
We see all the graphs are straight lines. Here we selected 
a = 0.9 because of the fact that small a in the weights of 
E 1 and E 2 is not a good estimate for d = 2 dimension. In 
Fig. 2 we have plotted ip against all dimensions obtained 
at the lowest pair and highest pair (-0 2 ) of length 
scales in our numerical calculation. The two broken-line 
graphs corresponding to the two ip calculated at two ex- 
treme length scales almost fall on the continuous graph 
showing a plot of 9 against same dimensions. 

Now we would like to make a few comment regarding 
the difference of our result (ip = 0) at the lower bound 
and that obtained in specially because we are u sing 
the simplified numerical scheme proposed in . In |ll| 
it has been tacitly taken that the barrier introduced or 
removed with the change in boundary condition passes 
through the nearest neighbour spins to the right hand 
corner one which is actually flipped to change the bound- 
ary condition. We argue that nobody knows where the 
domain wall will pass through corresponding to a change 
in the boundary condition and the domain wall width 
will actually be ignored if we concentrate on a sequence 
of spin flip as is done in Moreover, flip of a clus- 

ter of spins can be equivalently achieved by following se- 
quences of single spin flip but there is no guarantee that 
one would get through a barrier state of lower energy un- 
less a method of sequential excitation and stabilization 
is applied as has been explained in the introduction. All 
these things taken into consideration and specially when 
we do not know where the actual barrier will pass through 
or will be removed from when the boundary condition is 
changed, one must apply renormalization group approach 
to compare between the routes of going from one ground 
state to the other defined on appropriate length scales 
and the outcome is the same as that set by Fisher and 
Huse 0. 

I acknowledge discussions with J.K. Bhattacharjee, 
Jayanth Banavar and Amos Maritan. 
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